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ABSTRACT 

We present single-dish and interferometric millimeter line observations of the 
HH 111 outflow and its driving source. The physical conditions of the protostellar 
core have been determined from the emission of the millimeter line emission of CO 
and its isotopomers and CS with the IRAM 30m telescope, and the CO J = 7 — > 6 
line with the Caltech Submm Observatory. The molecular gas emission reveals a 
small condensation of cold (T = 20 — 25 K) and dense gas (n(H2) = 3 x 10 5 cm -3 ). 
The low-velocity outflowing gas has been mapped with the IRAM Plateau de Bure 
interferometer. The cold gas is distributed in a hollow cylinder surrounding the 
optical jet. The formation of this cavity and its kinematics are well accounted 
for in the frame of outflow gas entrainment by jet bow shocks. Evidence of 
gas acceleration is found along the cavity walls, correlated with the presence 
of optical bow shocks. The cavity has been expanding with a mean velocity of 
4 km s" 1 on a timescale of 8700 yr, similar to the dynamical age of the optical jet. 
The separation of the inner walls reaches 8" — 10", which matches the transverse 
size of the wings in the bow shock. CSO observations of the J = 7 — > 6 line show 
evidence of a high- velocity and hot gas component (T = 300 — 1000 K) with a low 
filling factor. This emission probably arises from shocked molecular gas in the jet. 
Observations of the 3 P 2 — 3 Pi [CI] line are consistent with C-type non-dissociative 
shocks. Mapping of the high-velocity molecular bullets B1-B3, located beyond 
the optical jet, with the IRAM PdBI reveals small structures of 3" x 7" flattened 
perpendicular to the flow direction. They are made of cold (T ~ 30 K), moderate 
density gas (n(H 2 ) = (0.5 — 1.0) x 10 4 cm -3 ). We find evidence that the bullets 
are expanding into the low-density surrounding medium. Their properties are 
consistent with their being shocked gas knots resulting from past time-variable 
ejections in the jet. 

Subject headings: ISM:jets and outflows - Stars:formation - ISMxlouds - 
ISM:individual:HH 111 



1. Introduction 

Herbig-Haro jets are manifestations of 
the accretion processes that are involved 
in the birth of stars, and as such represent 
a fossil record of the recent activity of their 
driving sources. The study of jets thus en- 
ables us to understand the past history of 
young stars as they build up their mass. As 
these jets penetrate their ambient medium, 
they transfer momentum and accelerate 
their surroundings and it appears that it 
is this interaction which causes the ubiqui- 
tous molecular outflows. For a recent re- 
view on these issues, see Reipurth & Bally 
(2001). 

One of the finest HH jets known is the 
HH 111 jet, located in the L1617 cloud 
in Orion (Reipurth 1989). It has a very 
highly collimated body, with a large num- 
ber of individual knots, and has large 
proper motions (300 to 600 km s -1 ). The 
jet flows at an angle of only 10° to the 
plane of the sky (Reipurth, Raga, & Heath- 
cote 1992). Beyond the collimated jet the 
flow continues, and it is in fact only a mi- 
nor part of a truly gigantic bipolar HH 
complex stretching over 7.7 pc (Reipurth, 
Bally, & Devine 1997a). Optical spectra 
are discussed by Morse et al. (1993a,b) 
and Noriega- Crespo et al. (1993). Near- 
infrared observations have revealed a re- 
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markable symmetry between the optical jet 
and a near-infrared counterjet, and that a 
second bipolar flow, HH 121, emerges from 
the source, suggesting that it is a binary 
(Gredel & Reipurth 1993, 1994; Davis, 
Mundt, & Eisloffel 1994, Coppin et al. 
1998, Reipurth et al. 2000). The optical 
jet is co-axial with a major well collimated 
molecular outflow (Reipurth & Olberg 
1991, Cernicharo & Reipurth 1996, Nagar 
et al. 1997, Lee et al. 2000). Beyond the 
visible jet, Cernicharo & Reipurth (1996) 
(henceforth CR) found a set of three high 
velocity CO bullets, which apparently rep- 
resent working surfaces that no longer pro- 
duce optical shock emission. The jet is 
driven by IRAS 05491+0247, a Class I 
source with a luminosity of about 25 L , 
surrounded by cold dust and gas (Reipurth 
et al. 1993, Stapelfeldt k Scoville 1993, 
Yang et al. 1997, Dent, Matthews, Ward- 
Thompson 1998). A small 3.6 cm radio 
continuum jet was detected around the 
source by Rodriguez & Reipurth (1994). 
Recently, infrared images obtained with 
HST have revealed that another source is 
located only 3 arcsec from the main one, 
which itself is a binary judging from the 
quadrupolar outflows; so altogether it ap- 
pears that the driving source is part of a 
hierarchical triple system (Reipurth et al. 
1999). In fact, the entire giant jet flow may 
be a result of the dynamical interactions 
between the three members of this triple 
system (Reipurth 2000). 

The short blue lobe found by Reipurth 
& Olberg (1991) appears to be a highly col- 
limated molecular jet about 90" long with 
a collimation factor of 9 (see their Figs. 2- 
3-4) in the 12" data of CR. The blue lobe 
ends at the bow shock P and the optical 
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knots Q,R,S. Nagar et al. (1997) showed 
an interferometric map at 6" resolution of 
the J = 1 — > line, which allowed them to 
resolve the structure of the low- velocity gas 
around the optical jet. The main features 
of the ionized and molecular gas in HH 111 
are summarized in Fig. 1. 

In this paper we present single dish 
and high resolution interferometric obser- 
vations of the body of HH 111, the pro- 
tostellar core and the high velocity CO 
bullets. We characterize the gas kinemat- 
ics and the physical conditions in the gas 
and report the detection of a hot tempera- 
ture component in the outflow. We discuss 
these results in relation to current theories 
on the coupling of HH jets and molecular 
outflows. 

We find that the molecular gas entrain- 
ment in the outflow is well accounted for in 
the frame of jet bow shock driven outflows, 
as modelled by Raga & Cabrit (1993) 
(henceforth RC93). Preliminary results 
were reported by Reipurth & Cernicharo 
(1995) and Cernicharo, Neri & Reipurth 
(1997). 

2. Observations 

2.1. Single-Dish Data 

The central source and the HH 111 out- 
flow were observed in the millimeter lines 
of CS and the CO isotopomers with the 
IRAM-30m telescope in 1990 and 1994. 
The CO observations have already been 
presented in CR. For these observations, 
we used as spectrometer an autocorrela- 
tor with 2048 channels and spectral res- 
olution of 37 kHz. The spectral resolution 
was degraded to tv 0.2 km s -1 . The re- 
sults are expressed in units of main-beam 



brightness temperature. The main-beam 
efficiency was measured to be 0.63 and 0.57 
at the frequency of the CS J = 2 — > 1 and 
J = 3 — > 2 transitions, respectively. 

The CO J = 7 -> 6 line at 806.651708 GHz 
was observed towards the central source 
and bullet B3 with the Caltech Submm 
Observatory (CSO) in December 2005. 
The receiver was tuned in double-side band 
and the system temperature ranged be- 
tween 5700 and 7500 K. This allowed us 
to observe simultaneously the 3 P2 — 3 Pi 
[CI] lines at 809341.97 GHz in the image 
band. An AOS was used as spectrometer, 
which provided a kinematical resolution of 
0.24 km s -1 . At this frequency, the beam 
size of the telescope is about 12" and the 
beam efficiency is 0.28. The calibration 
was checked towards Ori IRc2 in Orion A. 

2.2. Interferometric data 

The CO J = 1 — > emission along the 
HH 111 outflow was observed with the 
IRAM Plateau de Bure interferometer 
(PdBI). A mosaic of 5 fields was per- 
formed; the field centers were located at 
the offset positions (+10",0), (-10",+3"), (- 
30",+6"),(-50",+9") and (-70",+12") with 
respect to the source respectively : 05 ft- 51 m 46.2 
+02 o 48'29.5" (J2000). 0528+134 and 
3C120 were used as phase calibrators. An 
autocorrelator was used as spectrometer, 
with an element of resolution of 20 kHz, 
corresponding to a kinematic resolution of 
~ 0.06 km s^ 1 at the CO J = 1 -> fre- 
quency. The IRAM 30m data of the re- 
gion obtained by CR were used to provide 
the short spacings and combined with the 
PdBI visibilities. The combined map is 
shown in Fig. 5. The synthetic beam is 
3.79" x 2.43" at 2.6 mm and makes a P.A. 
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Fig. 1. — Map of the CO J = 1 — > outflow emission integrated between +0.25 km s -1 and 
+7.3 km s^ 1 obtained with the IRAM Plateau de Bure interferometer. We have indicated 
the main features detected in the molecular gas (the cavity and the high velocity bullets) 
gas and in the ionized gas (the optical jet, as observed with HST, the shocked knots and the 
large bow shock V). A dotted line traces the main axis of the jet. 



of —4° deg with respect to North. 

The CO J = 2 — * 1 emission of the 
high- velocity bullets Bl, B2 and B3 dis- 
covered outside the molecular cloud was 
observed with the IRAM PdBI in 1996. 
Three fields were observed centered re- 
spectively at 05 fe 51 m 34.65 s +02°48'50.50" 
(Bl), 05 h 51 m 31.12 s +02°48'57.50" (B2), 
05 A 51 m 28.06 s +02°49'08.02" (B3). The 
same phase calibrators as before were used. 
The resolution of the autocorrelator was 
set to 160 kHz, which provided a kine- 
matic resolution of ~ 0.16 km s -1 at the 
CO J = 2 — > 1 frequency. The synthetic 
beam is 1.96" x 1.26" at 1.3 mm and makes 
a P. A. of —192° deg with respect to North. 

To estimate the fraction of flux missed 
by the interferometer, we have compared 
the flux measured at the IRAM 30m tele- 
scope with the flux recovered by the in- 
terferometer towards bullet B2 in the ve- 
locity interval of maximum intensity in- 
tensity. At the IRAM 30m, the line in- 
tensity peak is 0.43 K, which corresponds 
to a total flux of 5.1 Jy, after correcting 
for the source convolution with the beam 
and taking into account the antenna ef- 
ficiency of lOJy/K at 1.3mm. At the 
same velocity (—36.8 km s -1 ), the flux in- 
tegrated over the reconstructed source was 
estimated with the "FLUX" procedure and 
found to be 4.95 Jy. We conclude that al- 
most all the flux has been recovered by the 



interferometer. 

3. The source 

Figure 2 shows the emission of the mil- 
limeter lines of CS and the CO isotopomers 
detected towards the VLA source and two 
nearby positions offset by 20" in right 
ascension. In the direction of the VLA 
source, double-peaked profiles are detected 
in all lines of various opacities, but the 
C 18 J = 1 -> and 13 CO J = 1 -> (see, 
e.g., CS J = 2 -> 1 and 13 CO J = 2 -> 1). 
In particular, the maximum of the singly 
peaked 13 CO J = 1 — > coincides with the 
8.4 km s _1 peak of the J = 2 — > 1 tran- 
sition. This allows us to exclude self- 
absorption as origin of the double peaked 
profiles detected. 

The 13 CO interferometric map of HH 111 
by Stapelfeldt & Scoville (1993) indicates 
a North-South velocity gradient in the di- 
rection of the core where the powering 
source of HH 111 has been formed. The 
CS map of Yang et al. (1997) taken with 
the Nobeyama interferometer also suggests 
such a velocity gradient. It has been in- 
terpreted by Yang et al. (1997) as due 
to infalling gas around the exciting source 
of HH 111. However, the fact that the 
emission from the CS lines show wings 
associated with the molecular outflow ren- 
ders difficult the interpretation of the data. 
Indeed, the gas kinematics is affected by 
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Fig. 2. — (bottom) Molecular emission of 
CO, 13 CO, C 18 CS and DCO+ towards 
selected positions along the optical jet of 
HH111, indicated above each upper panel. 
The dashed lines mark the velocity of the 
protostellar core (v = 8.4kms _1 ) and the 
cavity wall (v = 9.2 km s -1 ). 



several different physical processes, as ev- 
idenced by the presence of a well defined 
cavity (see below Sect. 4.3; also CR and 
Nagar et al. 1997). 

We do not detect any spatial shift in the 
velocity centroid of the molecular emission 
in the direction perpendicular to the out- 
flow. Hence, the double peaked line pro- 
files are not tracing the collapsing gas of 
the protostellar envelope, but most likely 
the interaction of the jet with the ambi- 
ent gas, constituting the rear part of an 
expanding cavity around the optical jet. 
This point is adressed more thoroughly in 
Sect. 4. 

The C 18 J = 1 -> transition traces 
the emission of the parental core at 8.4 km s _1 . 
The C 18 and 13 CO J = 1 -> were ob- 
served simultaneously by the same receiver 
(they lie in the same IF band), which of- 
fers the advantage of eliminating the un- 
certainties in the pointing and the rela- 
tive calibration. Towards the source, we 
measure main-beam brightness tempera- 
tures T 18 = 4.9 K and T 13 = 10.6 K for 
the C 18 and the 13 CO transitions, re- 
spectively. Both transitions are thermal- 
ized at the density of the core. Assum- 
ing a standard relative abundance ratio 
[ 13 CO]/[C 18 0] = 8, we derive an opacity 
t* 8 = 0.56 for the C 18 line. Analysis of 
the C 18 transition in the Large- Velocity 
Gradient approximation yields a gas col- 
umn density N(C 18 0) = 6.5 x 10 15 cm -2 
and a temperature of 15 K. Hence, the 
total gas column density of the core is 
N(H 2 ) = 3.3 x 10 22 cm' 2 . 

The CO J = 7 -> 6 to J = 2 -> 1 inten- 
sity line ratio measured at ambient veloc- 
ities, between +3.8 and + 11.5 km s -1 (see 
Sect. 4.1), is ~ 0.3, as expected for gas at 
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20 - 25 K and densities of 10 4 - 10 5 cm" 3 . 
We therefore adopt a value of 20 K in what 
follows. 

The dense gas tracers (CS) peak at 
vi sr = 8.4 km s" 1 in the central core re- 
gion. The lines are rather bright (T mb = 
3 K) and narrow (Av less than lkms" 1 ). 
Small wings are detected between 6.5 and 
llkms" 1 (Fig. 2). The emission draws 
a round condensation of 26" diameter 
(beam-deconvolved) or 0.06 pc. Analy- 
sis of the CS J = 3 -> 2 and J = 2 -> 1 
transitions in the Large- Velocity Gradient 
approximation yields typical densities of 
n(H2) = 3 x 10 5 cm" 3 and column densities 
N(CS) = 1.3 x 10 13 cm" 2 at the brightness 
peak. The lines are optically thin with 
r 21 = 0.44 and r 32 = 0.77. 

Emission in the 3 P2 — 3 Pi [CI] line was 
detected at the velocity of the ambient gas 
~ 9 km s _1 (Fig. 3; note that the [CI] line 
is in the image band). A gaussian fit to 
this emission yields an integrated flux of 
20 K km s" 1 , peak intensity of 1.5 K and a 
line width of 12kms _1 . Such a broad line 
width implies that the [CI] emission arises 
from the outflow and not only from the am- 
bient gas. This suggests a gas column den- 
sity of 1.0 x 10 18 cm" 2 for a temperature of 
20 K. The ratio of atomic carbon to car- 
bon monoxide column densities is therefore 
[CI] /CO ~ 0.3, a typical value measured in 
the Orion molecular cloud. 

4. The outflow 

4.1. Hot gas near the source 

No evidence of high-velocity wings was 
found towards the central source in the low 
J = l _, o and J = 2 -> 1 CO line profiles 
obtained with the IRAM telescope (CR; 
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Fig. 3. — (top). CO J = 7 — > 6 observed 
towards the central source (thick). We 
have superposed (thin) the J = 2 — > 1 line 
emission detected with the IRAM 30m 
telescope, (bottom) Variations of the 7- 
6/2-1 ratio R as a function of the ki- 
netic temperature and density. Calcula- 
tions were done for a CO column density 
of 1.0 x 10 17 cm" 2 in the LVG approxima- 
tion. The contours at (0.8,1.0,1.2) times 
the measured ratio R= 3.0 are drawn. 
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see also Fig. 3). CR pointed out that 
in case the outflowing gas would have a 
high excitation temperature, much above 
the 15-20 K estimated from the lowest two 
transitions, the opacity of the latter would 
be very low, making them extremely diffi- 
cult to detect. 

Hatchell et al. (1999) showed evidence 
for warm gas in the range 80 — 150 K in the 
HH 111 high- velocity bullets from observa- 
tions in the CO J = 4 — > 3 and J = 3 — > 2 
transitions. We have observed the central 
source in the high J = 7 — > 6 CO transi- 
tion at the CSO; the upper energy level lies 
at 155 K above the ground state, and is a 
much more sensitive probe of warm gas. 

The final spectrum is displayed in Fig. 3. 
The kinematical resolution was degraded 
to 2.4 km s _1 to improve the signal-to- 
noise ratio. The rms is 0.20 K per ve- 
locity channel. We have superposed the 
J = 2 — > 1 spectrum obtained with the 
IRAM 30m telescope, scaled to match the 
brightness peak of the J = 7 — > 6 line. The 
HPBW of both telescopes are similar at 
the concerned frequencies. We detect high- 
velocity wings both at blue and redshifted 
velocities, with a signahnoise of 4cr and 
above 5a, respectively. Emission is de- 
tected at velocities larger than 30 km s~ x 
with respect to the source. These wings 
are barely seen in the J = 2 — > 1 line. 
The J = 7 — > 6 line emission can there- 
fore be divided in 3 velocity intervals : the 
blueshifted wine; (-40 to +3.8 km s" 1 ), the 
"ambient" gas (+3.8 to +13kms _1 ), and 
the redshifted wing (+13 to +40kms~ 1 ). 

As illustrated in Fig. 3, the relative in- 
tensity of the CO J = 7 — > 6 with respect 
to the J = 2 — > 1 is much larger in the 
wings than in the ambient gas, which indi- 



cates highly different physical conditions, 
in particular a much higher temperature. 
We have estimated the physical conditions 
in the high-velocity wings from modelling 
both lines in the Large- Velocity Gradient 
approach. We have adopted the CO col- 
lisional coefficients determined by Flower 
(2001) for CO - ortho-H 2 collisions in the 
range 5 K-400 K. For temperatures beyond 
400 K, the collisional coefficients are ex- 
trapolated adopting a temperature depen- 
dence of ^(T/400 K). The signal : noise ra- 
tio in the wings is not high enough to study 
the excitation conditions as a function of 
velocity. Therefore, we consider the emis- 
sion integrated over the velocity interval 
assigned to each wing. We obtain a sim- 
ilar J = 7 — > 6/ J = 2 — > 1 line ratio ps 3. 
This value is probably a lower limit as the 
filling factor of the J = 7 — > 6 is probably 
less than that of the J = 2 -> 1. The LVG 
calculations show that, in order to account 
for a J = 7^6/J=2^1 line ratio of 
~ 3, the kinetic temperature has to be 
above 100 K at least, and the density larger 
than a few 10 3 cm~ 3 (see Fig. 3). In this 
range of density and column density, the 
J = 2 — > 1 transition is thermalized but 
the J = 7 — > 6 is not; the J = 2 — > 1 is op- 
tically thin while the J = 7 — > 6 is opti- 
cally thick. 

The velocity of this hot gas is very large, 
up to 150 km s -1 when deprojected. If the 
emission arises from entrained gas then the 
density of the hot gas is probably much less 
than in the cold low-velocity gas, where 
we estimate n(H 2 ) = 7 x 10 4 cm" 3 (see 
Sect. 4.2). The density of the hot CO gas 
is at most a few times 10 4 cm -3 , which im- 
plies a temperature of about 300 K up to 
1000 K or more. Spectroscopic diagnostics 
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in the optical range indicate typical densi- 
ties of 7000 cm -3 in the jet (Hartigan et al. 
1994). 

The CO column density of entrained gas 
in each wing lies in the range (0.3 — 1.0) x 
10 17 cm~ 2 . The mass of warm H 2 is about 
(0.3 — 1.0) x 10~ 2 M per beam, and the as- 
sociated momentum is 0.3 — 1.0 M & km s -1 , 
assuming a filling factor of 1. Integrated 
over the cavity, the mass of warm gas 
would amount to 0.02 — 0.07 M and the 
velocity-deprojected momentum would be 
~ 2— 7 M Q km s _1 . The overall momentum 
carried away in the hot CO gas would be a 
factor of 10 larger than in the low-velocity 
cold gas and the high- velocity molecular jet 
(Nagar et al. 1997; CR). We conclude that 
the filling factor of the hot component is 
probably much smaller than 1 . For a filling 
factor of 0.1, the same amount of momen- 
tum is carried away by the three compo- 
nents (molecular jet, high-velocity wings, 
low- velocity "cavity"). As a consequence, 
the hot CO emission is to arise from a nar- 
row region in/around the optical jet, and 
the size of the emitting region is a few arc- 
sec. However, we can not exclude that 
the emission arises from the whole cavity. 
If so, Kelvin-Helmholtz instabilities should 
grow along the jet. This has not been ob- 
served by HST. As discussed by Reipurth 
et al. (1997b; hereafter Rei97) the opti- 
cal jet propagates into a very-low density 
medium (n ~ 350 cm -3 ). Therefore, we 
favor the hypothesis that the hot CO gas 
emission comes from the jet itself. 

The spectra obtained by CR in the CO 
J = 2 — > 1 line further West in the flow, 
towards optical knots P to S, show simi- 
lar line wing profiles, typically 30 km s _1 
wide, with maximum intensity at the am- 



bient gas velocity, plus the signature of the 
molecular jet detected 5 — 8 km s 
wide isolated velocity component cen- 
tered at — 60kms _1 (see panels a 
and b in Fig. 3). Assuming similar ex- 
citation conditions in the jet, i.e. a 
J = 7^6/J=2^ 1 line intensity ratio 
of 3, we would expect a J = 7 — > 6 line 
intensity of ~ 0.4 K, which is below the 
detection limit allowed by our data. 

When averaging the J = 2 — > 1 line pro- 
files over the low-velocity outflow emis- 
sion region, between the source and off- 
set —70", a weak high-velocity component 
is detected at the 4<r level, which peaks 
at the same velocity as in the jet (see 
panel c in Fig. 4). Its narrow line width 
(Av ~ 3kms _1 ) is consistent with a mo- 
tion parallel to the optical jet (a spec- 
trum with 13 km s _1 kinematical resolution 
is displayed in Fig. 1 of CR). Assuming a 
minimum excitation temperature of 100 K, 
the mass of this component is 4 x 10~ 4 M & 
and the momentum carried away is about 
0.2 M & km s _1 , after correcting for velocity 
projection effects. 

High signahnoise J = 2 — > 1 data show 
that the velocity distribution between the 
molecular jet and the high- velocity wings is 
continuous (see e.g. the spectrum obtained 
near knot S in Fig. 4a). This supports the 
idea that the high-velocity wings consist 
of material entrained by/with the jet. ISO 
observations of molecular outflows have de- 
tected the emission of rotational lines of 
CO, H 2 and H 2 at temperatures between 
500 and 1200 K (Nisini et al. 2000; Lefloch 
et al. 2003). In the case of L1448-mm 
and HH 2, this hot gas emission appears 
to arise from extreme high- velocity clumps 
of 2" size along the molecular jet. 



8 



w 

III III u , , 


'"1 
i i i i 


I, 
I, 










: 'b') 1 1 1 


iii 

v l 
1 1 1 1 








:'cV ' 


1 1 1 1 

r ' 




I. 




i 



-100 -50 50 
V (km/s) 



Fig. 4. — High-velocity CO outflow emis- 
sion detected with the IRAM 30m tele- 
scope beyond the cavity, at 96" from the 
source (a), near optical knots P-S (b) and 
towards the cavity, within 70" from the 
source (c). 



As a conclusion, the hot CO high- 
velocity emission detected in HH 111 prob- 
ably traces gas that has been accelerated 
by internal shocks along the jet. 

4.2. The physical conditions 

The single-dish data allow to determine 
the physical conditions in the outflowing 
gas. The CO J = 2 -> 1 and J = 1 -> 
line intensity both peak at +9.7 km s" 1 
(Fig. 2). The gas at this velocity ex- 
tends along the optical jet. The emission 
maxima are at positions —10 to —20 arc- 
sec away from the exciting source of the 
HH 111 system; it appears as bright as the 
main core in the lines of CO and 13 CO. 
In the C 18 J = 1 — > line, a secondary 
component at +9 km s" 1 is detected from 
—30 arcsec to +60 arcsec from the exciting 
source along the outflow. A weak coun- 
terpart is detected at 7.7 km s -1 about 50" 
from the source, which coincides with the 
red outflow wing. The parameters of the 
secondary component were obtained from 
a gaussian fit to the C 18 lines, which 
yields a line width of 0.55 km s" 1 , a bright- 
ness temperature of 1.4 K, and an emission 
peak at +9.2 km s" 1 . This implies a col- 
umn density of ~ 2.5 x 10 21 cm" 2 provided 
the kinetic temperature is similar to that of 
the main core, as discussed above. We have 
estimated the density in the outflow wings 
from studying the secondary component in 
the CS lines. After removing the contribu- 
tion of the main outflow component to the 
spectra, we have determined from a gaus- 
sian fit main-beam peak temperatures of 
1.4 K and 0.75 K for the J = 2 -> 1 and 
J = 3 — > 2 lines from a gaussian fit. At 
a temperature of 20 — 25 K, this implies a 
gas density n(H 2 ) = 7x 10 4 cm" 3 , and a CS 
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Fig. 5. — Interferometric map of the CO J = 1 — > emission in HH 111. The optical HH jet, 
as observed with HST, is drawn in white. The synthetic beam (2.5" x 3.5" HPBW) is 
indicated in the lower-right panel. First contour and contour interval in the various panels : 
0.6 K in the first column; 0.6 and 0.9 K, respectively, in the second column; 0.6 K for the 
panels at 7.7 and 7.2kms _1 ; 0.6 and 1.0 for the panel at 6.4kms _1 ; 0.4 K for the panels at 
5.6 and 4.9 km s" 1 (third column); 0.4 K in the fourth column. 



column density N(CS) = 7 x 10 12 cm 2 . 

4.3. The PdBI observations 

We present in Fig. 5 the velocity channel 
maps of the CO J = 1 — > emission be- 
tween and +13kms _1 , obtained at the 
PdBI. The maps result from a mosaic of 
5 overlapping fields covering 120" x 40". 
Zero spacing data from the IRAM 30-m 
telescope has been merged with the PdBI 
data. The HPBW is 2.5" x 3.5" and is in- 
dicated in the lower-right panel. These ob- 
servations have a better angular resolution 
(ss a factor of 2), sensitivity and uv cover- 
age than the map of Nagar et al. (1997). 
We have superposed in white the emission 
of the optical jet, as observed with HST in 
Ha (Rei97). The data confirms the pres- 
ence of a hollow tubular structure with an 
outer diameter of 30" (2.1 x 10 17 cm), as 
previously proposed by CR and Nagar et 
al. (1997). Within 30" to the source, the 
geometry of the cavity is well fitted by 
a parabolic law. However, further away 
from the source, the gas distribution is bet- 
ter described by a hollow conical cylinder 
(Fig. 6) 

The front and rear sides of the cav- 
ity are detected between +2 to +8 km s _1 
and +8.5 to +12 km, respectively, hence 
at velocities much lower than the CO out- 
flow (— 60kms~ 1 , CR). Line widths in the 



blueshifted gas are rather large, with Av ~ 
4kms~ 1 at half-power. The walls of the 
cavity are seen particularly prominent at 
around +6kms _1 ; they present a wider 
separation as the distance from the source 
and the velocity increase (see Fig. 6). HST 
observations in the optical [SII] line indi- 
cate that the internal bow shock wings of 
the jet have a typical diameter of 5" (see 
Figs. 7-11 in Rei97) so that they fill only 
partly the molecular gas cavity. The walls 
extend as far as 60" away from the source 
(4.1 x 10 17 cm), beyond the HH jet. At 
—70", we measure an inner separation of 8" 
between the walls, which corresponds to an 
opening angle of 5 °. This separation is also 
the distance between the wings of the bow 
shock observed in the jet at this position. 
Emission at extreme velocities of +2.5 and 
+ 10.5 km s _1 is detected along the jet main 
axis, which provides a measurement of the 
expansion of the cavity around the jet. We 
derive an average radial expansion velocity 
of m 4kms~ 1 . This implies a kinematical 
age of 8700 yr for the outermost gas layers 
in the expanding cavity, at 15" away from 
the optical jet axis. This is similar to the 
kinematical age of the optical jet (CR). 

The velocity distribution on the rear 
side of the cavity provides strong con- 
straints on its geometry. The emission of 
the rear wall of the cavity is detected at 
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Fig. 6. — (top) Map of the CO J = 1 — > emission integrated between +0.25 and 
+7.3 km s -1 . The position of optical knots along the jet is marked with squares. The 
white line fits the cavity of gas swept-up in the propagation of bow shock V. Two dashed 
lines mark the opening of the cavity inner walls. The synthesized beam is drawn in the 
lower left corner. Contours range from 0.2 to 1.1 by step of 0.1 times 7.1 Jy/beam km/s. 
(bottom) Magnified view of the CO cavity. 



Fig. 7. — Velocity-position diagram of the CO J = 1 — > emission along cuts parallel to the 
HH jet. Positive (negative) cuts are made north (south) of the jet axis. Distance to the jet 
(in arcsec) is indicated in the upper left corner of each panel. First contour and contour 
interval is 0.5 Jy/beam (3.7 K). In the cut at —0.4" from the jet axis, we have indicated 
with thick white lines the front and rear cavity walls. In thick black dashed is marked the 
direction of the secondary CO component. 



about 9kms _1 , West of the source, where 
the outflow propagates at blueshifted veloc- 
ities (see panel at offset position (—20,0) 
in Fig. 2). The emission of the rear wall 
is characterized by relatively narrow line 
widths (1 — 2 km s" 1 ) and low- velocity mo- 
tions with respect to the cloud gas, much 
smaller than those detected on the front 
side. This implies that the rear side lies 
close to the plane of the sky, in agreement 
with the ~ 10 ° angle of the optical jet to 
the plane of the sky derived by Reipurth, 
Raga & Heathcote (1992) and the above 
opening angle of 5 °. The velocity distri- 
bution undergoes a positive gradient with 
increasing distance to the source, up to 
offset +50" along the jet. As the jet is 
blueshifted, this positive gradient is nat- 
urally explained if the tangential velocity 
is much larger than the component normal 
to the wall. 

We present in Fig. 7 the CO intensity 
distribution as a function of velocity in cuts 



made parallel to the jet direction. Emis- 
sion from the front cavity walls is detected 
from —11" to +11" from the jet axis. The 
location of the front and rear cavity walls 
is marked with thick white lines, in the 
cut made along the jet axis (offset —0.4"). 
The velocity shift between the front and 
the rear wall increases strongly with dis- 
tance x to the source, which implies that 
the tangential velocity also increases. The 
velocity field of the front and rear walls is 
well described by a linear law v oc x; the 
corresponding fits are displayed in Fig. 7. 

The maximum velocity (~ lkms -1 ) 
is observed ahead of the optical jet. A 
small gas knot is detected at the tip of 
the jet between offsets —45" (bow shock 
L) and —60" (knot O) in the panel at 
+1.4 km s -1 (see Fig. 5). A counterpart to 
this knot is detected in the redshifted gas, 
at +9.5 km s -1 . Emission at a lower ve- 
locity (+2.4 km s" 1 ) is detected upstream 
of the knot and off the jet axis (Fig. 5). 
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The whole wall appears to be accelerated 
along the optical jet up to this knot, as 
can be seen on the cuts at — 11" and +2.4" 
in Fig. 7). About 5" downstream of knot 
O the wings of bow shock P are detected 
(Fig. 8 in Rei97). Those wings are re- 
markably extended, 10", comparable to the 
width of the inner cavity "hole". The 
molecular gas has probably been acceler- 
ated in the wings of bow shock P. 

Near the source, the gas undergoes a 
strong acceleration, which terminates at 
the location of optical knot F, at 27" from 
the source (Figs. 4-5; Rei97). A symmet- 
rical feature displaying a strong gas accel- 
eration is also detected , which extends up 
to +15" East of the source. Both compo- 
nents are marked by a white dashed line in 
Fig. 7. 

We suggest that both features are trac- 
ing a recent ejection from the protostellar 
source, associated with knots E-F. This gas 
component extends from —6" up to +5" 
from the jet main axis (Fig. 7); hence the 
ejection takes place inside the large-scale 
cavity reported above. Indeed, based on 
STIS spectroscopic data of the jet obtained 
with HST, Raga et al. (2002) showed good 
evidence that optical knot F traces an in- 
ternal working surface resulting from time 
variability in the ejection velocity, on a 
scale of 250 yr. 

The velocity field appears much more 
complex (see panel at —0.4" in Fig. 7). It 
can be fitted by a linear law between the 
source and offset Aa = —20". One can 
note that the linewidth gets broader with 
decreasing distance to knot F, in panels at 
—3.2", +2.4" and along the jet. Between 
—25" and —40", the CO emission covers a 
broader range of velocities. Such broad- 



ening could trace gas accelerated by knots 
E-F, or in previous events like bow shock 
L. 

5. The CO bullets 

Beyond the optical bow shock V, CR 
discovered three CO bullets (hereafter re- 
ferred to as Bullet 1, 2 and 3, in order of 
increasing distance from the source) with 
separations of ~ 55". No optical emission 
is associated with these bullets. However, 
they are perfectly aligned along the direc- 
tion defined by the optical jet. Whereas 
the optical knots P and V are tracing 
shocked material, the molecular bullets 
Bl-2-3 have no heating source. Based on 
the association of optical knots with some 
high-velocity bullets, CR proposed that 
the latter bullets were tracing past shock 
episodes of HH 111. 

5.1. The PdBI observations 

High-angular resolution maps have been 
obtained with the PdBI interferometer in 
the CO J = 2^1 line (Fig. 8). The 
emission properties of all three bullets are 
very similar. The molecular emission from 
the bullets was marginally resolved by the 
IRAM 30m-telescope; it allowed to esti- 
mate a size of 8" — 10". The much bet- 
ter angular resolution of the PdBI observa- 
tions allows to resolve the emission. Each 
bullet is elongated in the North-South di- 
rection; the large-scale distribution is well 
described by an ellipse with typical size 
(beam-deconvolved) of 3" x 7" (2.0 x 4.8 x 
10 16 cm). We note that the transverse size 
of the bullets (~ 10") matches reason- 
ably well the inner diameter of the cavity 
around the jet. 



12 




Aa (") Aa (") Aa (") 



Fig. 8. — Interferometric maps of the velocity-integrated flux in the bullets Bl, B2, B3 in 
the CO J = 2 — > f transition. First contour interval are 0.10 and 0.04 Jy/beam respectively. 



Each bullet presents a complex kinemat- 
ical structure down to arcsec scale (Fig. 9). 
The line profiles are characterized by a 
broad plateau, spanning 10 — 15kms _1 , 
and a line core of 5 km s _1 HPFW. To 
study the gas kinematics in the bullets, 
we have considered separately the spa- 
tial distributions of the blue and red 
wings (Fig. 10). In all three bullets, the 
blueshifted component appears to lie ~ 1" 
(1 beam size) upstream of the redshifted 
component. It suggests a side- ways expan- 
sion of the bullets, which now propagate 
in a low-density medium, and are oriented 
at only a small angle out of the plane of 
the sky. 

5.2. The physical conditions 

Hatchell, Fuller & Ladd (1999) ob- 
served at the JCMT the CO J = 4 -> 3 
and J = 3 — > 2 emission near optical knot 
P, bow shock V and bullet B3. Under the 
hypothesis that the levels are populated 
at LTE and the lines are optically thin, 
they estimated a kinetic temperature of 



~ 100 K in the high- velocity ejecta. It is 
in bullet B3 that these authors derive the 
highest temperature, around 270 K. This 
is rather unexpected as B3 is the most dis- 
tant bullet from the source. We note that 
the transitions observed by the authors are 
not well suited to trace so high tempera- 
ture gas. Higher-J transitions of CO are 
more sensitive probes of the kinetic tem- 
perature. 

We have carried out an LVG analysis 
of the CO line fluxes measured towards 
B3, based on the data of CR and Hatchell 
et al. (1999). We have first corrected 
for the beam dilution, especially critical 
for the J = 1 — > transition, adopting a 
bullet size of 10", in agreement with CR 
and our PdBI observations (see below). 
We adopted an average line width Av = 
5.5 kms -1 , from the 30m and PdBI obser- 
vations. We conclude that it is not pos- 
sible to reproduce simultaneously all four 
transitions. The observed J = 4 — > 3 flux 
exceeds by a large factor (about 4) the in- 
tensity predicted by the model. It is in fact 
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Fig. 9. — Map of the J = 2 — » 1 emission in B3, integrated in velocity intervals of 0.4 km s 
Contour interval is 0.2 Jy/beam (1.87 K). 
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Fig. 10. — (bottom) Average spectrum of the bullets, in the region of intensity larger than 
50% the brightness peak, (top) blueshifted (redshifted) emission in solid (dashed) contour. 
First contour and contour interval are 0.10 and 0.04 Jy/beam, respectively. 



14 



Table 1: Physical parameters of the CO bullets mapped with the PdBI. 



Bullet 


Position 


Size 


Mass 


Density 




(J2000) j 


(") 


(1O" 4 M ) 


(10 4 cm- 3 ) 


Bl 


05:51:34.9, 02:48:51.1 


3.3 x 4.6 


3.1 


1.5 


B2 


05:51:31.4, 02:48:57.7 


3.3 x 6.7 


2.1 


0.59 


B3 


05:51:29.0, 02:49:08.0 


3.3 x 6.7 


1.9 


0.53 



the high flux measured in this transition 
which is responsible for the high temper- 
ature derived by Hatchell et al. (1999). 
However, there is no evidence for inter- 
nal shocks, which would heat the bullet 
material to such high temperatures. As 
discussed below, the cooling timescale of 
this gas is very short, so that its tem- 
perature should decrease much below the 
100 — 150 K measured in the optical knots 
P-V. We conclude that, most likely, the 
J = 4 — > 3 line intensity was highly over- 
estimated. As pointed out by Hatchell et 
al. (1999), the calibration of the 460 GHz 
receiver was highly uncertain at the time 
of the J = 4 — >■ 3 observations. 

If we repeat the analysis considering 
only the first three CO transitions and 
taking into account the uncertainty in the 
measured flux (about 10%), a good agree- 
ment is obtained for a CO column density 
N(CO) = 1.2 x 10 16 cm" 2 . The best agree- 
ment is then obtained for temperatures in 
the range 35 — 45 K and H 2 densities in the 
range 5 — 8 x 10 3 cm -3 (see below). 

Neufeld, Lepp & Melnick (1993,1995) 
have studied the radiative cooling of warm 
molecular gas T > 100 K in various as- 
trophysical environments. Applying their 
work to the case of the HH 111 bullets, we 
estimate a typical cooling time of 600 yr if 



the temperature of the bullets is as high as 
100 K. This is approximately the typical 
flight time between two subsequent molec- 
ular bullets (500 yr). The distance from 
bullet V, the last one for which we observe 
evidence of a shock, to bullet B3 is 120", 
that needs about 1500 yr to be crossed, al- 
most three times the CO cooling time. 

We conclude that, although the temper- 
ature of the bullets inside the molecular 
cloud is about 100 — 150 K, the molecu- 
lar bullets B1-B2-B3 are most likely made 
of colder gas, at a temperature of typi- 
cally 35 — 45 K. Under such physical con- 
ditions (ra(H 2 ) = 8000 cm" 3 , N(CO) = 
1.2 x 10 16 cm" 2 ), the J = 2 -> 1 transition 
is optically thin r ~ 0.3 and has an ex- 
citation temperature T ex ~ 20 K. This is 
slightly larger than the value derived by 
CR (15 K) although the other physical pa- 
rameters (mass, mass transport rate, mo- 
mentum) are essentially unchanged. 

From the above timescale, we estimate 
a time-averaged mass transport of 4 x 
10~ 7 M yr _1 , whereas Hartigan, Morse & 
Raymond (1995) derive a mass loss rate 
of 1.8 x lO _7 M yr _1 . The similar mass 
transport rates support the view that the 
CO bullets are not made of ambient mate- 
rial that has been accelerated to these very 
high velocities, but are mostly composed of 
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material ejected from the central engine. 
The momentum carried away by each bul- 
let amounts to 1.0 x 10 -4 M Q km s _1 yr -1 . 

Raga et al. (1990, 2002) proposed that 
optical knots form from time variable gas 
ejections in the jet. Such mechanism leads 
to the formation of "mini-bows" along the 
jet. In this framework, the high- velocity 
CO bullets would trace the previously 
shocked molecular counterpart of the opti- 
cal knots. It naturally explains the similar 
properties observed between the jet and 
the high-velocity bullets (size, kinematics, 
momentum). The weak high- velocity com- 
ponent detected between Bl and B2 (see 
Fig. 6 in CR) traces the molecular "inter- 
shock" gas ejected from the central source, 
between two "eruptions". 

6. Gas entrainment mechanism 
6.1. The cavity 

The kinematics of the outflowing gas is 
complex as it bears the signature of several 
ejection episodes in the optical jet, associ- 
ated with bow shocks at knots F and L 
(Rei97; Raga et al. 2002). We also find 
strong correlation between the low- and 
high- velocity outflow wings and the jet. 
It is therefore unlikely that a single fit of 
any outflow model is able to reproduce the 
kinematics of the gas observed along the 
whole flow length. As the eject a are lo- 
cated inside the cavity, close to the source 
(within ps 30"), we have decided to concen- 
trate on the global properties of the cavity, 
at larger distance from the source. Our 
interferometric data have revealed a large 
cylindrical cavity in which the g 
celerated following a "Hubble law" v oc x, 
distance along the jet axis. 



The formation of such large cavities is 
easily accounted for by the sweeping up of 
ambient material as a bow shock propa- 
gates into the parental cloud (RC93). In 
what follows, we show that that such gas 
entrainment models can indeed account for 
the molecular outflow emission in HH 111. 
We first determine the parameters of the 
"best" model that fits the HH 111 molec- 
ular gas distribution. 

As shown by RC93, the global shape 
of the cavity swept-up by an internal 
bow shock depends on two parameters : 
the density stratification of the ambient 
medium and the ratio between the dis- 
tance of the working surface to the source 
d and the standoff distance between the 
working surface and the bow shock apex 
r s . 

For a jet propagating into a medium 
of density p oc x' 3 , x measuring the dis- 
tance from the powering source along the 
jet axis, the following relation holds be- 
tween the jet cylindrical radius r c and x : 
r c oc x' 13 / 2 . As discussed above and in 
Sect. 4.3, the shape of the cavity near the 
source is suggestive of f3 = —1. 

The standoff distance r s between the 
ejection point in the working surface and 
the bow shock apex is related to the jet 
radius via : 



r, = r,- 



1/2 



where Q is the solid angle over which mass 
is ejected (Q ~ 27rv / 2), vq and Vj are the 
ejection and jet velocities, respectively (see 
Eq. 4 in RC93). 

In the case of strong working surfaces 
such as those observed in HH 111, the ra- 
tio between the jet and the ambient den- 
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sity pj/pa is ~ 100 (Raga & Binette 1991) 
and t>o — Vj in the adiabatic case. Hence 
r s /r-j ~ 6 x (vo/vj) < 6 and r s < 9", with 
a typical jet radius of 1.5" (Rei97). 

A simple determination comes from the 
shape of the cavity computed in the case 
(3 = —1. The cylindrical radius is typi- 
cally twice the standoff distance (see for 
example the top panel of Fig. 3); it is 
not very sensitive to the actual value of 
(3 nor the distance of the bow shock to the 
source. At zeroth order, the diameter of 
the cavity determined in the outer part is 

18" = 4r s , hence r s m 4.5". This simple 
determination is consistent with the upper 
limit obtained previously. It is also con- 
sistent with numerical simulations, which 
yield r s jr,j ~ 3. 

Since bow shock V lies at 148" from the 
source (Rei97), one has d/r s = 32. We dis- 
play in Fig. 6 the shape of the cavity swept- 
up by a bow shock as predicted for this set 
of parameters. This model fits nicely the 
walls of the cavity and the curvature of bow 
shock V. The agreement is very satisfying 
at Ao < —50". We have fitted the inner 
walls of the cavity by straight lines (drawn 
in dashed dark contours in Fig. 6); the op- 
tical jet emission appears to fill exactly the 
cavity walls up to Aa = —50". Extrapo- 
lating the fit up to bow shock V, we find 
that the separation of the walls at that po- 
sition is ~ 10". It matches well the size of 
the bow shock in the Ha optical line (see 
top panel in Fig. 6), as expected for freely 
expanding shocked material. 

For values of f3 of — 1, the model predicts 
that the velocity vs position along the axis 
is almost exactly linear over the whole flow 
length, except close to the source. This is 
again consistent with the Posit ion- Velocity 



diagram derived from our interferometric 
maps (see Fig. 7). 

The situation is more complex closer to 
the source : molecular gas emission is still 
detected beyond the fit, so that the agree- 
ment between the fit and the actual shape 
of the cavity is not very good. Several fac- 
tors may contribute to this. In particular, 
several ejections have taken place recently, 
which are now observed as knots E to L. 
Since this part of the cavity has experi- 
enced the crossing of two more bow shocks 
(P and L), the walls have been transferred 
additional momentum. This could con- 
tribute to our detection of entrained gas 
at further distances from the jet axis. 

These results differ somewhat from the 
previous conclusions of Nagar et al. (1997). 
Based on their interferometric map of the 
CO J = 1 — > line emission, these authors 
interpreted the velocity field and the spa- 
tial distribution of the molecular gas as ev- 
idence for an outflow entrained by a wide- 
angle radial wind, as modelled by Shu et 
al. (1991, 1995). The modelling of the 
Position- Velocity diagram of the emission 
along the jet is indeed in good agreement 
with the predictions of radial-wind driven 
outflows (Li & Shu 1996). As pointed out 
by Cabrit et al. (1997), a simple fit to the 
Position- Velocity distribution of the out- 
flowing gas does not allow to discriminate 
between jet-driven and wide-angle wind 
models, when the inclination with respect 
to the plane of the sky is low. 

We note that the physical conditions 
close to the source (less than 15") play an 
important role in constraining the param- 
eters of the wide-angle wind fitted by Na- 
gar et al. (1997). As can be seen on their 
Fig. 5, it appears that at larger distances 
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the radius of the cavity is constant rather 
than parabolic-like, as claimed by the au- 
thors, and the match between the observed 
velocity profile and the fit is less satisfying. 

In our opinion, the strong ejection ac- 
tivity near the source prevents conclusions 
about the role of a wide-angle wind in driv- 
ing; molecular g these ejecta are very 
efficient at transferring momentum to the 
ambient material. Interferometric obser- 
vations at a higher angular resolution than 
the data presented here are required to in- 
vestigate the relation the optical shocked 
knots hold to the molecular gas, and how 
the dynamics of the cavity is affected. 

However, the formation of the large- 
scale cavity observed up to bow shock P 
is well explained in the frame of the model 
of RC93. The jet entrains ambient molec- 
ular material swept-up in the propagation 
of bow shocks, which disturb the ambient 
gas at large distances from the jet axis. 
Our observations indicate that bow shock 
V played a major role in shaping the global 
cavity. 

6.2. The hot gas component 

Our observations of the CO 7-6 emis- 
sion bring evidence for hot outflowing gas 
in the cavity which is missed by obser- 
vations of the low-J transitions, due to 
the extremely low opacity of the latter. 
This is evidenced by the maximum radial 
velocity detected in the outflow wings in 
the J = 7 — > 6 and J = 2 — > 1 transitions 
(30kms _1 and lOkms -1 , respectively). It 
is therefore necessary to take into account 
the thermal structure of the outflow to dis- 
criminate between the various gas accelera- 
tion mechanisms that have been proposed. 

The hot gas component occupies prob- 



ably a much smaller region than the low- 
velocity cold gas, with a typical size of the 
emitting region of a few arcsec. The emis- 
sion likely arises from shocked gas in the 
jet, a conclusion supported by the high ve- 
locities and temperature measured in this 
gas. High signahnoise spectra obtained 
further downstream in the jet show a con- 
tinuous velocity distribution between the 
molecular jet and the high-velocity wings, 
and support a close association between 
both phenomena. [CI] observations do not 
show any evidence of CO dissociation close 
to the source; hence, the gas is probably 
accelerated in C-type, rather than J-type 
shocks, along the jet. 

Averaging the CO J = 2 — > 1 signal 
over the cavity reveals a faint molecular jet 
signature at GOkms" 1 , consistent with the 
velocity of the molecular features (bullets) 
detected further downstream (Fig. 4). The 
presence of such a hot component shows 
that there is no clear velocity separation 
between the outflowing gas and the jet it- 
self (see Fig. 4, and lower panel of Fig. 1 
in CR). Interferometric observations with 
the SMA in the CO J = 6 -> 5 transition 
would allow to determine the location of 
the emission. Observations of the high-J 
CO lines with the Herschel Space Obser- 
vatory would allow to clarify the physical 
conditions in the outflowing gas. 

7. Conclusions 

We have presented single-dish and inter- 
ferometric line observations of the HH 111 
outflow and its driving source. The physi- 
cal conditions of the protostellar core were 
determined from the emission of the mil- 
limeter lines of CO and its isotopomers, 
the millimeter lines of CS, all of them ob- 
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served with the IRAM 30m telescope, and 
the CO J = 7 — > 6 line observed with the 
CSO. The molecular gas emission reveals a 
small condensation of cold (T = 20 — 25 K) 
and dense gas (n(H 2 ) = 3 x 10 5 cm -3 ). 

The low- velocity outflowing gas has 
been mapped with the IRAM Plateau de 
Bure interferometer. It is distributed in a 
hollow cylinder that surrounds the optical 
jet. The formation of this cavity and its 
kinematics are well accounted for in the 
frame of the model of RC93. The jet en- 
trains ambient molecular material swept- 
up in the propagation of bow shocks, which 
disturb the ambient gas at large distances 
from the jet axis. Our observations in- 
dicate that bow shock V played a major 
role in shaping the global cavity. The gas 
velocity is found to vary almost linearly 
with distance to the source over the whole 
flow length, except close to the source. The 
cavity has been expanding with a mean ve- 
locity of 4 km s _1 on a timescale of 8700 yr, 
similar to the dynamical age of the opti- 
cal jet. The separation of the inner walls 
reaches 8 — 10" near knots P and S, which 
corresponds to the transverse size of the 
wings in these bow shocks. 

The more recent ejections detected close 
to the source are still efficient at transfer- 
ring momentum to the ambient material. 
Regions of strong gas acceleration are de- 
tected along the jet, which are correlated 
with the optical bow shocks at knots F and 
L. 

Observations of the high-excitation line 
CO J = 7 — > 6 reveal high-velocity out- 
flowing gas (v = 150 km s _1 after deprojec- 
tion) near the central source. This high- 
velocity component is hot (T = 300 — 
1000 K) and fills only a small fraction of 



the beam. We propose that this emission 
arises from shocked gas in the optical jet. 

Emission of the molecular bullets Bl-2- 
3 was mapped with the IRAM PdBI in the 
CO J = 2 — > 1 transition. The emission is 
distributed in structures of 3" x 7" flattened 
perpendicularly to the jet main axis. Their 
kinematics is complex but an overall side- 
ways expanding motion is detected inside 
all of them. Their physical properties have 
been reassessed; we find a low tempera- 
ture for the bullets (~ 30 K), in agreement 
with the cooling timescale for shocked gas 
initially heated at temperatures of 100 K 
or beyond, and cooling down while fly- 
ing freely. Their density lies now in the 
range (0.5 — 1.0) x 10 4 cm -3 . These molec- 
ular bullets probably form in the shocks re- 
sulting from time-variable ejections in the 
powering source, in the manner proposed 
by Raga et al. (2002) to account for the 
presence of mini bow shocks and knots in 
the optical jet. In such a scenario, the 
weak high-velocity molecular jet detected 
between the molecular bullets represents 
intershocked gas ejected from the central 
source. 
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